Transition Temperature of a Magnetic Semiconductor with Angular Momentum j 
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We employ dynamical mean-field theory to identify the materials properties that optimize T c for 
a generalized double-exchange (DE) model. We reach the surprising conclusion that T c achieves a 
maximum when the band angular momentum j equals 3/2 and when the masses in the rrij = ±1/2 
and ±3/2 sub-bands are equal. However, we also find that T c is significantly reduced as the ratio of 
the masses decreases from one. Consequently, the search for dilute magnetic semiconductors (DMS) 
materials with high T c should proceed on two fronts. In semiconductors with p bands, such as the 
currently studied Mn-doped Ge and GaAs semiconductors, T c may be optimized by tuning the band 
masses through strain engineering or artificial nanostructures. On the other hand, semiconductors 
with s or d bands with nearly equal effective masses might prove to have higher T c 's than p-band 
materials with disparate effective masses. 



The discovery of dilute-magnetic semiconductors 
(DMS) with current transition temperatures above 170 
K[l-3] initiated an active search for the optimal mate- 
rial for spintronic device applications [4]. However, the 
ferromagnetic transition temperature of magnetic semi- 
conductors involves many parameters that are difficult to 
control experimentally. To facilitate the search for new 
magnetic semiconductors with high transition tempera- 
tures, we evaluate T c for a double-exchange system with 
general angular momentum j using dynamical mean-field 
theory (DMFT). Surprisingly, T c is found to reach a max- 
imum for j = 3/2 and for equal light and heavy hole 
masses. However, we also find that T c decreases signifi- 
cantly when the ratio of the masses is reduced. 

The notion of using magnetic semiconductors in spin- 
tronic devices dates back to the 1960's, when europium 
chalcogenides [5] and chromium spinels [6] were exten- 
sively studied. Before the appearance of the latest gen- 
eration of III-V DMS grown by molecular beam epitaxy 
techniques [1, 2], II- VI [7] and IV- VI [8] DMS were de- 
veloped by alloying non-magnetic semiconductors with 
magnetic ions. Improvement in growth techniques has 
pushed the T c of Gai_ K Mn x As to values above 170 K 
[3]. As we discuss later, the large value of T c in delta- 
doped GaAs [9] might be associated with the reduced 
magnetic frustration when the Mn ions are restricted to 
two-dimensional planes. 

In our calculation we employ the DMFT, which was 
formulated in the late 1980's by Miiller-Hartmann [10] 
and Metzner and Vollhardt [11]. It has since developed 
into one of the most powerful many-body techniques for 
studying electronic models such as the Hubbard [12, 13] 
and DE [14-19] models. Since DMFT becomes exact in 
the dilute limit, it is a good starting point for studying 
DMS. Recent work on DMS materials has used DMFT 
to study variants of the DE model [20, 21] with less than 
one local moment per site. A DE model with one local 
moment per site and large coupling constant J c provides 
an upper limit to the transition temperature for a sys- 
tem with exchange coupling between the local moments 



and charge carriers. Perhaps more importantly, the be- 
havior of n/j < x holes in the impurity band of a DMS 
with x < 1 Mn atoms per site is very similar to that of a 
DE model with filling p — rih/x < 1 [21]. So as long as 
J c is sufficiently large to produce a well-defined impurity 
band, the qualitative results of this model should be in- 
dependent of the precise magnitude of J c /W [19]. As we 
shall see, a generalized DE model with one local moment 
per site and large J c also has the distinct advantage that 
analytic results are possible for any angular momentum 
j of the charge carriers. 

In semiconductors like GaAs, the angular momentum 
of the holes is obtained from the vector sum of the s = 
1/2 spin with the / = 1 orbital angular momentum of 
the p bands. The j = 3/2 band lies highest in energy 
while the spin-orbit split j = 1/2 band lies an energy 
A so ps 340 meV below [22]. Consequently, almost all 
of the holes in Mn-doped GaAs populate the j = 3/2 
band, which in turn is split by crystal fields [23] into a 
rrij = ±3/2 sub-band with heavy holes and a rrij = ±1/2 
sub-band with light holes. More exotic semiconductors 
with (I = 2) d bands, such as the chalcogenides [5], may 
contain carriers with total angular momentum j = 5/2. 

The optimization of T c is performed in two stages. 
First, we calculate T c for a DE model where the charge 
carriers have angular momentum j but equal masses in 
all 2j + 1 sub-bands. After concluding that the system 
with j = 3/2 and a half-filled lower or upper band has 
the highest T c , we examine the effect of different band 
masses in the rrij = ±1/2 and ±3/2 sub-bands. 

The Hamiltonian of a generalized DE model with carri- 
ers (holes or electrons) of angular momentum j and equal 
masses in all sub-bands is given by 

k i,k,k' 

(1) 

where cj^ and Ck Q are the creation and destruction op- 
erators for an electron with angular-momentum compo- 
nent rrij = a (a — —j, — j + 1, . . . , j) and momentum k, 
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Si = Srxii is the spin of the local moment (treated clas- 
sically) at site Rj, and J Q ^ are the (2j + l)-dimensional 
angular-momentum matrices. Repeated spin indices are 
summed. 

Within DMFT, the local effective action A e g(m) is 
governed by the bare Green's function Go, which con- 
tains dynamical information about the hopping of elec- 
trons onto and off neighboring sites. Because A c s(m) 
is quadratic in the field variables, the full local Green's 
function G{w n )af3 may be readily solved by integration, 

with the result [14] G(iv n ) = ((G (w„) _1 + J c J-m) 1 ) m , 
where v n = (2n + 1)ttT are Matsubara frequencies and 
underlined quantities are matrices in (2j + 1) x (2j + 1) 
spin space; (X(m)) m = J dfl m X(m.)P(m) denotes an 
average over the orientations m of the local moment with 
a probability -P(m), and J c = J C S. This relation is solved 
for a semicircular density-of-states with full bandwidth 
W in terms of the electron filling p (p — 1 corresponds 
to one electron per site so that < p < 2j + 1) and 
the local- moment order parameter M = (rn z ) m , which 
becomes nonzero below T c . 

As J c increases, the (2j + l)-degenerate band splits 
into 2j + 1 sub-bands, each labeled by quantum number 
rrij and centered at energy —nijj c . Due to the effect of 
electronic correlations, the full bandwidth of each sub- 
band is lowered from W to W = W/y/2j + 1. So prior 
to taking the limit of large J c for the rrij sub-band, we 
must rewrite the chemical potential as \x — ~rrijj c + Sfi 
where \Sfj,\ < W'/2. 

Since matrix quantities like G(w„), G (iv n ), and the 
self-energy £(w„) = G (w„) _1 — G_(iis n )~ 1 have al- 
ready been averaged over the orientations of m, they 
are spin diagonal and may be expanded in powers of J_ z 
as A = A I_ + A\J_ Z + A 2 J^ + . . . where A n is of or- 
der M n . So to linear order in M, only the first two 
terms in this expansion are retained and the bare inverse 
Green's function may be parametrized as Gjj{iv n )~ x = 
(z n - nij J c + R n )l + Q n J_ z where z n = iv n + Sfi. 

For large J c and to linear order in M, we solved self- 
consistent equations for R n and Q n . The Curie temper- 
ature satisfies the condition 



while the filling p = T ^ TrRe\G(ii/ n )] 
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These two expressions generalize previous results [14, 18] 
for j = 1/2. Quite naturally, a system without partially 
filled sub-bands (p an integer) has a vanishing Curie tem- 
perature because the carriers are unable to hop to neigh- 
boring sites without incurring an infinite cost in coupling 
energy. 

For a given j, the largest T c always occurs in the sub- 
bands with rrij = ±j because those holes or electrons are 
able to most effectively take advantage of the exchange 
coupling that mediates the ferromagnetism between the 
local moments. For j = 1/2, 3/2, 5/2, 7/2, and 9/2, the 
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FIG. 1: Ferromagnetic transition temperature for a double- 
exchange model with general angular momentum j. (a) The 
dependence of T c /W on filling for j = 1/2 (solid black), 3/2 
(dotted red), 5/2 (dashed purple), 7/2 (long-dashed orange), 
and 9/2 (dotted green), (b) The maximum transition tem- 
perature T™* x /W (in a half-filled band with rrij = ±j) versus 
(2j — l)/2 with the limiting value l/(kny/2j + 1) plotted in 
the dashed curve. Large j values, although unphysical, are 
plotted to illustrate the asymptotic behavior. 



dependence of T c on filling within the rrij = ±j sub-bands 
is plotted in Fig. 1(a). As expected, T c is particle-hole 
symmetric within each sub-band and is the same for a 
system with p electrons (1 — p holes) or p holes (1 — p 
electrons) per site. So the largest Curie temperature, 
T c max , is obtained when the lowest or highest sub-band is 
half-filled. 

For a half- filled sub-band with rrij = ±j, we obtain 
an analytic expression for T'™ ax by converting the Mat- 
subara sum into an integral (assuming that T c max /VF is 
small), with the result 
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This expression gives a slightly larger value than the 
exact result obtained from Eq.(2), which is used to plot 
T c max /VF versus (2j - l)/2 in Fig. 1(b). For large j, 
T max _ (W'/4tt){1 - 7r/4j} saturates at W'/in, as in- 
dicated by the dashed line in Fig. 1(b). Compared to 
the maximum T c max for j = 1/2 of 0.0219VF, T c max for 
j = 3/2 of 0.0256VF is 16% higher. Even systems with 
j = 5/2 and 7/2 have higher T c 's than one with j = 1/2. 
The suppression of T c for j = 1/2 is due to the stronger 
fluctuations at the smallest angular momentum. 

To interpret Fig. 1(b), keep in mind that in the ab- 
sence of magnetic impurities, a semiconductor is charac- 
terized by the bandwidth W of the conduction band and 
by the angular momentum j of the charge carriers. So for 
a class of undoped materials with the same bandwidth, 
the highest transition temperature will be found in the 
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FIG. 2: Cartoon illustrating the magnetic frustration due to 
the chirality of the holes mediating the ferromagnetism. (a) 
The minimum energy configuration for a pair of Mn moments, 
(b) the frustration that occurs when 4 Mn atoms lie at the 
corners of a tetrahedron. 



magnetic semiconductor with j — 3/2. A weakness of 
the present approach is that the local moments with spin 
S are treated classically whereas the charge carriers with 
angular momentum j are treated quantum-mechanically. 
Including the fluctuations of the local moments will fur- 
ther suppress T c . Therefore, the optimum DMS system 
will have S j so that fluctuations of the local moment 
can be neglected. This assumption is barely satisfied in 
Mn-doped GaAs, where S = 5/2 [3]. 

Having established that the semiconductor with j = 
3/2 has the highest Curie temperature, we now examine 
the effect of different masses in the heavy mj = ±3/2 and 
light rrij = ±1/2 sub-bands of the parent material. For 
GaAs, the band masses are nih — 0.50m and mi — 0.07m 
(m is the electron mass) with a ratio r — mi/rrih = 0.14 
[22]. As pointed out by Zarand and Janko [24], the differ- 
ent band masses introduce magnetic frustration because 

the kinetic energy K = J2w, a p ^.a/sc^a ^ is onl y di " 
agonalized when the angular momentum of the charge 
carriers is quantized along the momentum direction (for 
r = 1, the quantization axis is arbitrary). Due to the 
chirality of the holes mediating the ferromagnetism, the 
moments of a pair of Mn atoms at sites Ri and R2 pre- 
fer to align perpendicular to the unit vector ri2 along the 
R2 — Ri direction, as sketched in Fig. 2(a). 

This anisotropy is easy to understand. Most of the 
carriers that couple Mn atoms 1 and 2 will have mo- 
mentum k along ri2 and angular momentum j parallel 
to the Mn moments. If the Mn moments lie parallel to 
1*12, only heavy holes with angular-momentum compo- 
nent j • ri2 =3/2 would couple the two Mn. But if 
the Mn moments lie perpendicular to ri2 as in Fig. 2(a), 
holes with angular momentum perpendicular to r±2 but 
momentum parallel to T12 would contain all four angular- 
momentum components j -r^ = ±3/2 (heavy) and ±1/2 
(light). Due to the higher mobility of the light holes, 
the latter coupling is more effective. When four or more 
Mn atoms are distributed in space, the Mn interactions 
will be frustrated because some of the pairs cannot obtain 
their lowest energy configuration. This is shown schemat- 
ically in Fig. 2(b), where the circles (crosses) indicate 
pairs of moments that are (not) in their minimum en- 
ergy configurations. The large value of T c in delta-doped 




FIG. 3: Results for a double-exchange model with j = 3/2 
and different carrier masses in the heavy and light sub-bands, 
(a) The interacting density-of-states for j — 3/2 and a range 
of ratios r between the effective masses; at the right hand 
corner the non-interacting densities-of-states for r = 0.2 and 
0.6. (b) The dependence of T c /W t on filling for j = 3/2 and 
the same values of r. 



GaAs [9] may be related to the reduced magnetic frustra- 
tion when Mn ions restricted to two-dimensional planes 
can order with their moments perpendicular to the plane. 

The calculation of T c is performed using two semi- 
circular densities-of-states such as the ones sketched in 
the right hand corner of Fig. 3(a) with bandwidths Wi 
and Wh — rWi. As the interaction J c is turned on, the 
interacting density-of-states separates into four identi- 
cal sub-bands which are not particle-hole symmetric. In 
Fig. 3(a), we plot one of these sub-bands for a range of 
effective mass ratios r. For r < 1, the interacting density- 
of-states is always weighted towards lower energies, where 
the heavy holes dominate. 

In Fig. 3(b), we plot T c /Wi versus filling for the same 
values of r used in Fig. 3(a). Not surprisingly, T c de- 
creases as r decreases from I due to the magnetic frus- 
tration introduced by the chirality of the charge carriers. 
Hence, the maximum T c for a fixed light bandwidth is 
obtained when the band masses are identical. 

As r drops from I to 0, T™ ax /Wi plotted in Fig. 4 
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FIG. 4: Relation between T c max and the ratio of effective 
masses r = mi/rrih- In the main panel, T c max /W; is plot- 
ted versus r = mi/irih, with the experimental ratios denoted 
for several semiconductors. The horizontal arrows denote 
r c mB /W| for j = 1/2 and j = 5/2 (equal masses in all sub- 
bands) . The inset displays the filling for which the largest Tc 
is obtained versus mi /nth- 

decreases by 50% from 0.0256 to 0.0127. For the ex- 
perimental value r = 0.14 [22], T c max w 0.0133I4 7 ; has 
decreased 48% from its unfrustrated value. The experi- 
mental parameters r for a variety of semiconductors with 
j = 3/2 are also indicated in Fig. 4. Based on numer- 
ical results with Mn concentration x < 1 [21], however, 
the suppression of T c may not be quite so dramatic as 
suggested by this figure. Notice also that a system with 
j = 1/2 or j = 5/2 and nearly identical band masses 
may have a higher transition temperature than one with 
j = 3/2 but very different effective carrier masses. 

In the inset to Fig. 4, we plot the filling at which T c is 
maximized versus the ratio of masses. This filling varies 



from p = 0.5 at r = 1 to roughly 0.74 at r = 0. So for 
r < 1, the transition temperature is no longer particle- 
hole symmetric within the lowest sub-band. This may ex- 
plain why the transition temperature in Mn-dopcd GaAs 
continues to increase for hole fillings supposedly larger 
than half of the Mn doping [25, 26]. Based on this work, 
we expect that T c reaches a maximum for a filling close 
to 70%. 

To summarize, we have examined the general depen- 
dence of the transition temperature of a magnetic semi- 
conductor on the angular momentum and filling of the 
charge carrier bands. For a fixed bandwidth of the par- 
ent compound, T c is maximized when j = 3/2, the same 
angular momentum carried by the holes in GaAs and Gc. 
For j > 3/2, the angular momentum of the charge carri- 
ers effectively suppresses T c due to the narrowing of the 
impurity band by electronic correlations. However, the 
suppression of T c due to the different effective masses in 
p-band semiconductors is even more significant. Because 
of the strong reduction of T c by the magnetic frustra- 
tion, s-band (j = 1/2) or c?-band (j = 5/2) semiconduc- 
tors with more nearly equal band-masses may prove to 
have higher T c 's than p-band (j = 3/2) semiconductors 
with disparate effective masses. For the class of j = 3/2 
semiconductors, T c will be maximized when the effects 
of magnetic frustration are reduced as much as possible. 
This may be achieved either by applying strain, thereby 
altering the band structure, or through digital doping [9]. 
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